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The non-perturbative electron-positron pair production (Schwinger effect) is considered for space- 
and time-dependent electric fields E(x,t). Based on the Dirac-Heisenberg- Wigner (DHW) formal- 
ism we derive a system of partial differential equations of infinite order for the sixteen irreducible 
components of the Wigner function. In the limit of spatially homogeneous fields the Vlasov equation 
of quantum kinetic theory (QKT) is rediscovered. It is shown that the quantum kinetic formalism 
can be exactly solved in the case of a constant electric field E(t) = Eq and the Sauter-type electric 
field E(t) = Eq sech 2 (t/r). These analytic solutions translate into corresponding expressions within 
the DHW formalism and allow to discuss the effect of higher derivatives. We observe that spatial 
field variations typically exert a strong influence on the components of the Wigner function for large 
momenta or for late times. 

PACS numbers: 12.20.Ds, 11.15.Tk 



I. INTRODUCTION 

Pair production in strong external electric fields is 
in many respects aparadigmatic phenomenon in quan- 
tum field theory [l|-y|. It is nonperturbative in the 
coupling times the external field strength. It exempli- 
fies the nontrivial properties of the quantum vacuum, 
as it manifests the instability of the vacuum against 
the formation of many-body states. In general, it de- 
pends strongly on the spacetime structure of the exter- 
nal field, such that the pair production process is ex- 
pected to exhibit features of nonlocality, final state cor- 
relations and real-time dynamics. Moreover, it is a non- 
equilibrium process in quantum field theory and as such 
belonging perhaps to the least-well understood branch 
of modern field theory. Whereas pair proliferation is ex- 
pected to occur at the critical Schwinger field strength 
E c = m 2 /e ~ 1.3 x 10 18 V/m, recent studies have sug- 
gested that pair production might become observable al- 
ready at lower but dynamically modulated field strengths 
|3-Ho|. These estimates of the required field strengths 
indicate that pair production might already become ac- 
cessible at future high-intensity l aser sy stems such as the 
extrem e lig ht h 
XFEL [Mill. 

Computing pair production in a complicated space- 
and time-dependent field such as a high-intensity pulse is 
by no means straightforward. Many different theoretical 
methods, such as the propertime method [3J , WKB tech- 
niques |16l - l2G | , the Schrodinger- Functional approach [21| , 
functional techniques |22|.|23|. quantum kinetic equations 
24 - J29fl , being also closely related to scattering techniques 
30( , various instanton techniques I31M35B , Borel summa- 
tion [36( , propagator constructions |37| , and worldline 
numerics |38| have been developed to study pair produc- 
tion in external fields. Most of those approaches have 
only been applied to one-dimensional temporal or spa- 
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tial inhomogeneities, see (39| for the only true multidi- 
mensional case. Also, finite-temperature contributions 
have been determined which under the assumption of lo- 
cal thermal equilibrium first occur at the two-loop level 
[4(J. l4l|. For thermal pair production from more general 
initial states, see |42h44| . 

For both, a profound understanding of the phe- 
nomenon as well as reliable quantitative predictions for 
realistic cases, a formalism that can deal with arbitrary 
space- and time-dependent fields is urgently required. 
This is also stressed by recent observations of charac- 
teristic and potentially easy to detect signatures of pair 
production in the momentum distribution of the pairs 
which has turned out surprisingly sensitive to the sub- 
cycle structure of high- intensity pulses (43, |46| , also ex- 
hibiting information about the quantum statistics of the 
particles involved [43, [48|. Such a formalism based on 
suitable real-time correlation functions is indeed avail- 
able and has already been studied in the context of pair 
production [42 . The present work is devoted to explor- 
ing this DHW formalism, putting it into the context also 
of other work such as quantum kinetic equations, and 
performing first systematic studies with the aid of both 
exactly soluble cases and within approximative schemes. 

This paper is organized as follows: In Sec. [IT] we briefly 
review the Dirac-Heisenberg- Wigner (DHW) formalism, 
adopting already a notation which will prove to be ad- 
vantageous in the following. We describe how the Quan- 
tum Kinetic Theory (QKT) enrerges as a specific limit 
of the DHW formalism and present some analytical so- 
lutions. In Sec. Mil we introduce a derivative expansion 
and discuss its region of validity. In Sec. IIVI we con- 
clude and provide an outlook. Details about the QKT 
are summarized in App. |XJ The analytical results for 
the irreducible components of the Wigner function in the 
constant electric field and Sauter-type electric field are 
given in App. [B] 



II. THE EQUAL-TIME DHW FORMALISM 

A classical statistical one-particle system is described 
by probability distributions J-{x, p; t) in 6-dimensional 
phase space {x, p}. The generalization for a relativis- 
tic quantum field theory is obtained by choosing an 
appropriate density operator and performing a Wigner 
transformation to 8-dimensional phase space {x^,^}. 
The corresponding Wigner operator W(x,p) is manifestly 
Lorentz covariant but the associated Wigner function 
(n,\W(x,p)\£l) may not have a clear physical interpre- 
tation [5 (J . l51| . 

Alternatively, one may drop the manifest Lorentz co- 
variance in favor of a canonical time evolution from 
the beginning and start with an equal-time density op- 
erator. The corresponding Wigner operator W(x,p;t) 
is then defined in 6-dimensional phase space {x,p}. 
It is an advantage of this approach that the Wigner 
function (il\W(x, p; t)\fl) might be interpreted as quasi- 
probability distribution in analogy to classical physics. It 
is an additional benefit that the equation of motion might 
be formulated as initial value problem |49l . l52fl . Alterna- 
tively, one could also start with the Lorentz covariant 
formulation and switch to the equal-time formulation by 
performing an energy integral over po [53H55J . 

We will adopt the equal-time formulation throughout 
this paper. Due to the fact that we have dropped mani- 
fest Lorentz covariance anyway, we will also fix the gauge 
from the beginning. We will choose the temporal gauge 
Aq = throughout, such that the electric and magnetic 
fields are calculable from the vector potential A(x, t) ac- 
cording to: 

E(x, t) = -d t A{x, t) , B(x, t) = V X A(x, t) . (1) 



will allow for introducing a properly defined kinetic mo- 
mentum variable p. In terms of the center-of-mass coor- 
dinates x = ^(^l + X2) and s — xi — X2 , it reads: 



C(x;s;t) = e ICJ -U2 



/_'' ; 2 A(x+\S,t)-sd\ 



[tt(a 



(3) 
The Wigner operator is then defined as the Fourier trans- 
form of C(x; s;t) with respect to the relative coordinate 
s, such that the arguments are the center-of-mass coor- 
dinate x, the kinetic momentum variable p and time t: 



W(x,p; t) = -±fd 3 s e-^ 3 C{x; s; t) 



(4) 



Note that if we had defined Eq. ([2} with ^ (x~2 , t) instead 
of \l/(a^2,t), the corresponding Wigner operator would 
have been Hermitian. With our definition, W(f,p, t) is 
not Hermitian but transforms like a Dirac matrix: 



VV t (x,p;i)=7°>V(£,p;i)7 



(5) 



In general, the Wigner function is then defined as the ex- 
pectation value of the Wigner operator (il\W(x,p;t)\Q.) 
with respect to the full interacting vacuum. However, due 
to the fact that we are mainly interested in describing 
Schwinger pair production in the following, we restrict 
ourselves to the vacuum state in the Heisenberg picture 
|O) = |0): 

W(x,p;t) = -±fd 3 se- i P-s(0\C{x';s]t)\0) . (6) 

In order to derive the equation of motion for the Wigner 
function, we take the time derivative of Eq. © and take 
the properly gauge fixed Dirac equation 



A. Derivation of the DHW formalism 

In this section we define the equal-time Wigner opera- 
tor W(x, p; t) in the presence of an external electromag- 
netic field. By applying a Hartree approximation for the 
electromagnetic field, i.e. treating it as C-number field 
instead of a operator- valued quantum field, we are able 
to derive the equation of motion for the corresponding 
Wigner function (£l\VV (x , p; t)\tt). 

For this, we consider the following equal-time density 
operator of two Dirac field operators in the Heisenberg 
picture: 

C(fi,^;t)^e- 8e A lAV '^ d5 '[vI/(xl,t),*(x 2 ,i)] , (2) 

where we have dropped the Lorentz indices for simplic- 
ity. Here we choose the equal-time commutator, since the 
equal-time anticommutator is trivially fulfilled for spinor 
fields. Additionally, in order to preserve gauge invariance 
we include a Wilson line factor with an integral of the vec- 
tor potential over a straight line. In fact, the choice of 
the integration path is not unique, but the present choice 



ij°d t + ry • V* - ieA(x, t) - m) *(x , t) = (7) 



into account. In the course of the derivation we adopt 
a Hartree approximation of the electromagnetic field, 
which should be a good approximation for high field 
strengths. This means that we replace the operator- 
valued electromagnetic quantum field by a C-number 
electromagnetic field: 

(0|i^(£,i)C(x;s;<)|0) — > F^(x,t)(0\C(x; s; t)\0) . (8) 

Diagrammatically, this approximation corresponds to ig- 
noring higher- loop radiative corrections. Physically, this 
implies that final state interactions as well as mass shift 
effects are ignored. This derivation finally yields the 
equation of motion for the Wigner function: 

D t W = -\f) s [ 7 °7, W] - im [7 , W] - iP { 7 °7, W} , 
with D t , Dg and P denoting the following non-local 



pseudo-differential operators: 

,1/2 
A = d t +e d\E(x + iXVf;,t)-Vjt , 

7-1/2 

r 1 / 2 
A = V s +e d\B(x + iXVft,t) x V ? , (10) 

.7-1/2 

/- 1 / 2 ., 

P = p - ie d\\B(x + i\Vp,t) xV ? . 

7-1/2 

As the Wigner function W(x,p;t) is in fact a Dirac 
matrix, we may expand it in terms of irreducible compo- 
nents by choosing an appropriate complete basis set of 
4x4 matrices {1, 75, 7 M , 7^75,(7^}. Actually we choose 
16 real functions (from now on called DHW functions) 
which transform under orthochronous Lorentz transfor- 
mations as scalar &(x,p; t), pseudoscalar p(x,p; t), vector 
Un(x,p;t), axialvector a /i (x,p;i) and tensor ty, v (x, p\ t) , 
respectively: 



j»{x,t) 



/"'J 



W(x, p; t) = - [Is + i7 5 p + 7 /i v /i + 7 M 75a M + ct^ v \ 

(11) 
Inserting this decomposition into the equation of motion, 
Eq. @ , and comparing the coefficients of the basis ma- 
trices, we find a PDE system for the 16 DHW functions. 
Introducing the compact notation for the tensorial com- 
ponents: 

(tfl)^ = toi - tjo , (i 2 ) i = Cijktjk , (12) 

this system reads: 



As 
Ap 

A wo + A? • v 

A &o + A? ' a 

Aw + A? ^o 

A a + A? ao 



-2P-ti = (13) 

+ 2P-t 2 = 2m a (14) 

(15) 

= 2m p (16) 

+ 2P x a = -2m t*i (17) 

+ 2P x v = (18) 



A *i + A x t 2 +2Ps 
Ati 2 -Ax ti - 2P P 



2m v (19) 

(20) 



Note that for spatially homogeneous electromagnetic 
fields F^ v {x : t) = F^ v {t), an enormous simplification oc- 
curs as the non-local operators Eq. (fTU|) reduce to local 
ones: 



A = d t + eE(t) ■ Vp 
A = V^ + eB(t) x V- 
P = p 



(21) 



It has been shown previously [49j| , that some of 
the DHW functions can be given an intuitive inter- 
pretation, whereas others do not have a classical ana- 
logue. First, the symmetrized electromagnetic current 



[\l/(x, t),7 M \l/(a?, £)] |0) is expressed as 
d 3 p 



ffrt) 



(2ny 



■^(x,p;t). 



(22) 



Additionally, several conservation laws concerning phys- 
ically observable quantities like the total charge Q, the 
total energy £ , the total linear momentum V and the 
total angular momentum M. are valid, 



j t {Q^,V;M} 







with: 



Q = e / dT v (x, p; i) , 

£ / dT \p ■ v(x,p]t) + ms(x,p;t)] 

\E(x,t)\ 2 + \B(x,t)f 
V = I dTp^j (x,p;t) 

d 3 xE(x,t) x B(x,t), 



d 3 x 



(23) 



(24) 



(25) 



(26) 



M 



dT 



xx pv (x,p;t) - -si(x,p;t) 



+ d 3 xxx E(x,t)x B(x,t) , (27) 

with dT = d 3 xd 3 p/(2ir) 3 denoting the phase-space vol- 
ume element. According to these expressions, we may 
associate s(x,p;t) with a mass density, vo(x,p;t) with 
a charge density, w(x,p]t) with a current density and 
a(x, p; t) with a spin density. Another important con- 
servation law concerns the norm of the Wigner function 
itself: 



- J dT Tr [W(x , p; t) W f (x, p; t)] = , (28) 

which translates into a conservation law for the 16 DHW 
functions. 



B. Quantum kinetic theory (QKT) as limit of the 
DHW formalism 



In this subsection we show that the DHW formalism 
in the case of a spatially homogeneous, time-dependent 
electric field E(x,t) = E(t)e^ and vanishing magnetic 
field B(x, t) = yields the well-known Vlasov equation of 
QKT for Schwinger pair production |24l. l2a. l27j| . For this, 
we first calculate the Wigner function for pure vacuum 
to obtain appropriate initial conditions. In a second step 
we simplify the PDE system (03) - ((20j) to an ODE sys- 
tem [43, which turns out to be equivalent to the Vlasov 
equation [56f . For an analysis of the relation between the 



Wigner function and QKT for several examples of pair 
production in non-abelian fields, see |57| . 

In order to calculate the Wigner function for pure vac- 
uum W vac (x,p; t), we consider first the general expression 
Eq. (J5]) for vanishing vector potential: A(x,t) = 0. We 
first decompose the Dirac field operator in its Fourier 
basis 



*(#,*) 



d 3 q 

(2tt) 3 



mty*- 



(29) 



and introduce a decomposition in terms of anti- 
commuting creation/annihilation operators as well as 
four-spinors 

Hq, t)=J^u s (q, t)a s (q) + v 3 (-q, t)b\ (-q) . (30) 



Evaluating the vacuum expectation value and taking ad- 
vantage of the four-spinor completeness relations, we fi- 
nally obtain for the vacuum Wigner function: 



W vac (x,p; t) = - 2^p7 t lm _ ^ ■ $ 



(31) 



with u>(p) = y/m 2 +P 2 . Comparing this expression with 
Eq. (fTTj) . we immediately see that (a) in the pure vac- 
uum only 4 DHW functions do not vanish and (b) these 
vacuum functions do not depend on x and t: 



S>vac \P) 
^vac® 



1m 
2p 



(32) 
(33) 



After fixing the vacuum initial conditions, we consider 
next the PDE system Eq. CE3]) - $ZU\) for E(x, t) = E(t)e 3 
and B(x, t) = in more detail: Due to spatial homogene- 
ity, the DHW functions do not depend on the variable x 
and hence all spatial derivatives vanish. As an immediate 
consequence, wq (p; t) decouples completely. Additionally, 
due to the fact that the DHW functions {p, ao, t2}(p; t) 
are subject to a closed set of equations which does not 
couple to the non- vanishing vacuum initial conditions, 
these functions have to vanish as well. As a consequence, 
the PDE system for former 16 DHW functions reduces 
to a PDE system for the remaining 10 DHW functions 
w(p;t) = (s,v,a,ti)(p;i): 



[8t + eE(t)d P3 ] w(p; t) = M(p)w(p; t) 



(34) 



Here, w(p; t) is a column vector and A4(p) is the following 
10 x 10 matrix: 



M(f) 



itli 



/ ° 


2p T 





-2p x -2m 





-2p x 


\-2p 


2m 


( 


~P3 P2 \ 


r =\ 


P3 -pi 


\ 


-Pi P\ / 



(35) 



(36) 



The PDE system Eq. (|3~4"|) will be simplified by apply- 
ing the method of characteristics. We introduce a new 
parameter a and assume that the originally independent 
variables depend on this new parameter: 



p = Tx[a) and t = r{a) 



(37) 



Imposing the following equality for any function J-{p; t) 
depending on the former independent variables p and t: 



jL + eE(t)J- HP\ t) = -LH^), r(a)) 



(38) 



we find a = t = t and n(q,t) = q— eA(t)e%. Note that 
7?((f, t) denotes the time-dependent kinetic momentum on 
a trajectory, whereas <f, which serves as an integration 
constant in the method of characteristics, corresponds 
to the canonical momentum. Additionally, we still have 
the notion of a phase-space kinetic momentum p. These 
three types of momenta have to be clearly distinguished 
in the following. To be consistent throughout this paper, 
we always denote: 

p kinetic momentum in phase space , 

q canonical momentum , 

7? (<f, t) kinetic momentum on a trajectory . 

On the one hand, any function defined in phase space pos- 
sesses only an explicit time dependence and will hence- 
forth be denoted by F(p]t). On the other hand, func- 
tions depending on the time-dependent kinetic momen- 
tum 7? (if, t) show both an explicit and an implicit time 
dependence and will be denoted by J-(q,t). 

Formally, the method of characteristics is applied to 
the PDE system Eq. (f3~4"]) by replacing p by ir(q,t), such 
that the relation between the phase-space DHW func- 
tions and the DHW functions on a trajectory reads: 



w(<f,i) = ™(p\t)\ p _ > .f_ e x w 
w{p;t) = w(q,t)\ q ^ f+eA{t) 



(39) 

(40) 



Consequently, the PDE system Eq. (|3~4"|) becomes an 
ODE system, with the former time-independent matrix 
M (p) becoming a time-dependent quantity A4 (<f, t) : 



— w(q,t) = M(q,t)w(q,t) . 



(41) 



In order to proceed, we seek an appropriate basis to span 
w(q,t), such that Eq. (|4Tj) reduces to a simple form: 



10 



t(q,t) = -2Y,X i {q,-t%(q,t) , 



(42) 



with the factor —2 chosen for later convenience. To this 
end, we exploit the vacuum initial conditions Eq. (1321) 
(|3"3")l and choose the first basis vector ei (q, t) such that 
in pure vacuum the first coefficient xiac(<f;^vac) = 1, 



whereas all other coefficients Xvac(9^vac) vanish. Con- 
sequently, we find a subset of basis vectors: 



<Bi(q,t) 



<B 2 (q,t) = 



1 



e±u(q,t) 



<B3(<f,*) 



/ m \ 

u(q,t) 

V o J 



mir 3 (q 3 ,t^ 

n(q,t)n 3 (q 3 ,t) -w 2 (q,t)e 3 



d 



(43) 



f_ 



7r(q, i) x e 3 
-me 3 / 



with e ± = ^/m 2 + ^ and u(q, t) = i/ e i + 7r 3(?3,i), 
which form an orthonormalized, complete set: 



M(q,t){ e 2 

©3 

(Pi 



«3 
-CE2 



As a consequence, only the coefficients X^ 1 ' 2 '^*?^) 
couple to the initial vacuum state whereas all other co- 
efficients X l {?it) vanish. This means that w(q,t) is fully 
characterized by: 



){q,t) = -2 y £,X i {q,t)ii(q,t) 



(46) 



Next we introduce f(q,t) = 1 — X 1 ®*) parametrizing 
the deviation from the vacuum state, such that in pure 
vacuum f v &c(q,t va _ c ) = 0. Additionally, we define: 



UJ 2 (q,t) 



(47) 



If we consider the ODE system Eq. (|4"Tj) together with 
the relations Eq. ([44]) - (|4"5"]l. we obtain: 



dt 



Mt) = Q(<f,*) x 2 (?,i) , 



(48) 



dt 



dt 



X 2 {q,t) = Q(q,t)[l- f(q,t)]-2u(q,t)x 3 (q,t) , (49) 



X 3 (?,t) = 2£(<f,i)x 2 (9,i) 



(50) 



together with vacuum initial conditions / va c(<?, ivac) = 
Xvac(9,ivac) = Xvac(<f^vac) = °- Tnis 0DE system is 
nothing but the well-known Vlasov equation of QKT in 
its differential form [56| (cf. also App. \K§. Note that 
f(q, t) denotes the single-particle momentum distribution 
function in quantum kinetic theory. Thus, the DHW 
formalism in the presence of a spatially homogeneous, 
time-dependent electric field E(x,t) = E(t)e 3 is com- 
pletely equivalent to QKT. However, the DHW formalism 
is the more general approach since it allows for any time- 
and space-dependent electromagnetic field whereas the 
Vlasov equation is restricted to spatially homogeneous, 
time-dependent electric fields. 

For some special cases, an exact solution of QKT can 
be found (see Sec. IHC|) . such that we are able to calcu- 
late the DHW functions as well. To this end one uses the 
representation Eq. (j46|) and projects back to phase-space 
Eq. (|40p . For Schwinger pair production in spatially ho- 
mogeneous, time-dependent electric fields, for instance, 
one finds that only 7 of the possible 16 DHW functions 
contribute: 



s(P;t) 

v±(p;t) 

v 3 (p,t) 

ai (p; t) 

a 2 (p;t) 

ti,3(P;t) 



2m , ,_, _ s 2m p 3 



w(p) 



xHfrt) 



e±w(p) 



x\m 



-—rrx l {p;t) 
u(p) 



±P3 2/-,n 

x (pit) 



ej_uj(p) 



2p 3 



u(p) 



x l (p;t) 



2e ± 2 



u(p) 



x'ipit) 



2p 2 . 3 



e_L 



X%Plt) 



2 Pi ,,3 



e± 
2m 



x J (p;*) 

x 3 (p;t) 



(51) 
(52) 
(53) 
(54) 
(55) 
(56) 



C. Exactly solvable electric fields 

In this subsection we derive the analytic expres- 
sions for the single-particle momentum distribution func- 
tion f(q,t) of QKT for both the constant electric field 
E(t) — E and the Sauter-type electric field E(t) = 
Eq sech (t/r) . The construction of the solution can be 
oriented along the lines of QKT, cf. App. O We first 
seek an analytic solution for g^ (q, t) of Eq. (|A6[) and 
determine its normalization such that it coincides at 
t vac -> -oo with G (+) (q, t) defined in Eq. (|A17|) . Accord- 
ing to Eq. (|A22|) . we are then able to calculate the Bo- 
goliubov coefficient /3(q,t) and, consequently, the single- 
particle momentum distribution function f(q,t) accord- 
ing to Eq. |A20)) . 

As soon as this solution is know, we are able to calcu- 
late the non- vanishing coefficients X 2 (<?; t) accord- 
ing to Eq. (|4"B"| - (f5TJ|) . As an immediate consequence of 
Eq. ([ST]) - (|56"|) , all the non- vanishing DHW functions can 
be calculated as well. 



1. Constant electric field 

A constant electric field E(i) = Eq might be repre- 
sented by the vector potential: 



■Eat 



(57) 



A{t) 
such that Eq. (|A6[) reads: 

i + («s + eS t) 2 + ze^o) 5(«s,*) = • (58) 



(<l 



Note that for notational simplicity we do not explicitly 
indicate the dependence on the orthogonal canonical mo- 



mentum via e\ 



qj_. Dynamically, q^ is the only 



relevant paramater such that the situation becomes ef- 
fectively one-dimensional. Introducing the dimension- 
less parameter r\ = e^/eEn and performing the variable 
transformation: 



q 3 + eE a t 



we see that g(u) will only depend on one dimensionless 
variable u due to the linear relation between q^ and t. As 
a consequence, the differential equation Eq. (|58p turns 
into the parabolic cylinder differential equation ([58|. 
Chapter 19): 



[dl + \v? + \{i + V )\~g(u)=Q 



(60) 



This second order differential equation has two standard 
solutions which are given by: 



g(+\u) = N^D_ 1+iTl/2 (-ue- 
5<->(u) = N^D_ iv/2 (-ue^ ) 



(61) 
(62) 



with D v (z) being the parabolic cylinder function and 
N^' being normalization factors. For u — > — oo, these 
solutions behave asymptotically as: 



g {+) {u) u 

g('\u) 



-oo i 



- + 2 1 og (|„|) + f ] 



II— 7— DO 



-i[v+iio g (|«i; 



(63) 
(64) 



On the other hand, the adiabatic mode functions are 
given by, cf. Eq. (|A17|) : 



G {±) (u) = 



^i&(u ,u) 



i 2 e v / 2i] + u 2 (yj2ii + u 2 T u) 



(65) 



with: 



e = E /E c 



(66) 



The dynamical phase Q(uq,u) can be explicitly calcu- 
lated as soon as we fix uq. Here, we choose the symmetric 
point uq = 0, such that the definite integral, 



0(0, u) 



1 



(59) yields: 



du'y/2q- 



(67) 



\j2i-j + u 2 + 2r\ log (u+ \J2r\ + u 2 \ - i] log 2ij 



(68) 



Consequently, we fix the normalization constants N^ > 
such that g^(u — > -oo) = G^^u — > — oo): 



AT(+) = 
7V(-) = 



1 



p z[j[l+log(2/r,)]- 



4j e T 



\fi- 



-il[l + io & {2/ n )] e -^ 



(69) 

(70) 



According to Eq. (|A20[) and Eq. (|A22|) . we are then able 
to calculate the single-particle momentum distribution 
function. Taking into account the general relation for 
parabolic cylinder functions: 



d z D u {z) = \zD v (z)-D 1+v (z) 
we finally obtain: 



(71) 



/(«) 



V 2 ')+« 2 



2-q + u 2 -u) D_ 1+iv / 2 (-ue 1 1 ) - 2e l * D ir] / 2 {-ue l 2) 



(72) 



Finally, we may show that we obtain the Schwinger result 
for the pair production rate in a constant electric field if 
we consider the limit u — > oo. To this end, we take the 
leading term in the asymptotic expansion of the parabolic 
cylinder functions in Eq. (|72|) . Neglecting terms of the 



order 0{u 1 ), they are given by: 



D-l- ir ,/2(— u**) 

D_ tri/2 (~u l2 i) 



-* r(&) e ' <+3 g '" e s 



,(73) 
(74) 



such that the asymptotic behavior of f(u) is given by: 
lim ]{u) = 2e _7n7 . (75) 

u— >oo 

As a consequence, the Schwinger pair production rate 
per volume and time n[e + e~], i.e. the first term in the 
Schwinger expression for the vacuum decay probability 
Sill, is found: 



[e + e-]=lAd t Mt) = ^ 



(76) 



which completes our analytical solution for the constant 
electric field. 



2. Sauter-type electric field 



The Sauter-type electric field E(t) 



might be represented by the vector potential: 

A{t) = -E t tanh (1) , 
such that Eq. (|A6j reads: 



E sech (t/r) 

(77) 



d? 



(q 3 + eE Tt&nh(±)) + 

+ieE sech 2 (^)}g(q 3 ,t) = 0.(78) 



Again, we only indicate the dependence on q 3 whereas 
the dependence on the orthogonal canonical momentum 
will not be denoted explicitly. In the following, we intro- 
duce the dimensionless variable: 



= i[l + tanh(i)] 



(79) 



such that t — > — oo corresponds to u — > whereas t — > oo 
corresponds to u — > 1. Additionally, we introduce the 
Keldysh parameter 7 = l/(me7") and dimensionless mo- 
mentum variables q 3 = q 3 /m and k 2 = q^/m 2 , such that 
the dimensionless kinetic momentum on the trajectory 
^3(93? u ) an d the dimensionless energy variable uj(q 3 ,u) 
read: 



^3 (93, u) 



_ q : ,+eEgT(2u-l) _ 



= <73 



2u-l 



w 2 (g 3 ,w) = l + K 2 +7r|(g 3 ,u) 



(80) 

(81) 



Within these new variables, the differential equation 
Eq. J78]) reads: 



4 7 2 e 2 u(l - u)d u {u{\ - u)}d u + Q 2 {q 3 , u) + 

+4ieu{l-u)]g(q 3 ,u)=0 . (82) 



In order to solve this differential equation, we apply an 
ansatz for g(q 3 ,u): 



g(g3,u) 



= W3_0) 



3(432°)' 



*r* {1-u) 1 *r« h(q 3 ,u), 



(83) 



Plugging this ansatz into Eq. (|82l) yields the hyper- 
geometric differential equation ([58|, Chapter 15) for 
h(q 3 ,u): 



u(l — u)d 2 + (c — [a + b + l]u)d u — ab h(q 3 , u) = 



with 



(84) 



a(Q3) 
6(93) 



(* 



S(q 3 ,0) _ g(q 3 ,l) 
ye V 7 ' 2 2 



l + ^i-^ + ^M, 



c(ffi) = 1-£D(&,0). 



Note that these parameters do not depend on u. The two 
linearly independent solutions h^(q 3 ,u) in the neigh- 
borhood of the singular point u = are given by: 

h^{q 3 ,u) = A (+) F(5,&,5;w), (86) 

xF(l-a,l-b,2-c;u) ,(87) 



with -F(a, 6, c; u) denoting the Gauss hypergeometric 
function. Taking Eq. (|83p into account, the asymptotic 
behavior of g^'(q 3 , u) for u — > is: 



(+),„ v u->0 + -i "'' 3 '"' log(«) 

k '(q 3 ,u) — > e 2 ~<' &y ' 



9 



g { } (q3,m) 



»-»o + el £^io g («) 



(89) 



On the other hand, we are again able to give an analytic 
expression for the adiabatic mode functions: 



G^\q 3 ,u) = 



= =F«9(<33,«0,«) 



v /2m 2 uj(q 3 , u)[Q(q 3 , u) T n 3 (q 3 , u)] ' 



(90) 

Again, the dynamical phase ®(q 3 ,iio,u), which is given 
by the following integral: 



©(<Z3, U ,U) 



U , , -y/l + K 2 +S§(lj3,M') 



7e J tin 



27c Ju 



u' (1 — u') 



(91) 



can be analytically calculated. For uq 7^ {0, 1} and u — > 
0, this phase splits into a relevant divergent part and an 
irrelevant regular part ty(q 3 , ito, 0): 



&(q 3 ,u ,0) = y(q 3 ,u ,0) 



3(g3,0) 
27c 



log(u) , (92) 



such that the normalization constants N^ is fixed ac- 
cording to g^\q 3 ,0) = G^(q 3 ,0): 



N V 



o =Fi*(«3."0.°) 



^2m 2 2(<j3,0)[2(<j3,0)T5f3(<33,0)] 



(93) 



Again, the single-particle momentum distribution func- 
tion is calculated according to Eq. (|A20p and Eq. (|A22|) : 



/(<73, u) = N f (q 3 ) (l + g|^l) | [Q{q 3 , u) - (1 - «)w(«3, 0) - uu(g 3 , 1)] F(5, 6, 5; u) - 2i 7 e U (l - u)d u F(a, b, 5; «) 



(94) 



r 



with the normalization factor Nf(q 3 ) being given by: 



Nf(q 3 ) 



2O(g 3 ,0)[S(g 3 ,0)-ff 3 (g3,0)J 



(95) 



and 



<9 U F(5, b, c; u) = ^F(l + a, 1 + 6, 1 + c; z) . (96) 



Similarly to the constant electric field, we may give a 
simple expression for the asymptotic single-particle mo- 
mentum distribution function. Applying a linear trans- 
formation formula for hypergeomctric functions and con- 
sidering the asymptotic limit u — V 1, one first obtains: 



/(53,1) 



"(g3,l)+'*3(g34) 



D(<j3, 0)D(g3,l) S(g3,l)— #3(93,1) 



a& r(i+c)r(i+a+&-c) 
£ r(i+5)r(i+6) 



.(97) 



Applying the transformation formulae for gamma func- 
tions 



T(l + a) = aT(a) and |r(l + i&)| 2 



sinh(7r6) ' 



(98) 



we obtain a relatively simple analytic expression for the 
asymptotic single-particle momentum distribution func- 
tion /((?3,1): 

2smh(^[^+Q(q 3 A)-Q(q 3 ,0)])sinh(^ 7 [^-Q(q 3 ,l)+Q(q 3 ,0)]) 



sinh(^S(g 3 ,l)) sinh(^£((j3,0)) 



(99) 



III. 



INFLUENCE OF A SMALL SPATIAL 
INHOMOGENEITY 



In this section, we discuss the influence of a small 
spatial inhomogeneity along the direction of the time- 
dependent electric field E(x, t) based on the analytic re- 
sults for both the constant electric field and the Sauter- 
type electric field. In order to estimate the effect of higher 
derivatives, we adopt a derivative expansion and deter- 
mine the ratio between the first derivative and higher 
derivatives. Note that we again ignore the effect of mag- 
netic fields for simplicity, B(x,t) = 0. 



A. Derivative expansion 

We consider a space- and time-dependent electric field: 
E(x>,t) = E(t)[l + A(x 3 )}e 3 , (100) 



where |A(x3)| <C 1 describes a small deviation from the 
spatially homogeneous electric field. The equation of mo- 
tion for the 16 DHW functions w(x,p;t) read: 

D t w(x,p;t) = M(Vg,p)w(x,p;t) , (101) 



with wf(x,p;t) and Ai(Vs,p) satisfying Eq. (|13[) 
(p?0")) . We perform a derivative expansion of the pseudo- 
differential operator D t in Eq. (JTUJ), such that: 



D t = d t +eE(t)d P3 +eE(t) 



A(x 3 )d p 



A"(x 3 ) a3 
24 P3 



(102) 

where A" (0:3) denotes the second derivative with respect 
to x 3 . We may consider w(x, p; t) in an expansion as well: 



w(x,p;t) = w (0) (p;t) + w ( - 1) {x,p;t) 



(103) 



with w(°\p;t) being the exact result for the case of a 
spatially homogeneous electric field and w^ 1 ' (x,p;t) de- 
noting a small deviation from the zeroth-order solution 
due to the spatial inhomogeneity. Neglecting terms of 
the order A(x 3 )w^' (x,p;t), we obtain: 



d t + eE(t)d P3 - M (V s ,p)\ w« (x, p; t) « 

-eE(t) \A(x 3 )d P3 - ^P±d 3 J w^(p;t) , (104) 



such that the spatially homogeneous solution w^ (p; t) 
acts as a source term for the spatially inhomogeneous 
solution. In order to estimate the parameter regime for 
which the omission of derivatives higher than linear might 
be a good approximation, we consider a simple model: 



A(x 3 ) = 5 cos(f) , 



(105) 



with L being the length scale of the spatial variation and 
Sq <C 1 being its amplitude. Introducing the dimension- 
less variable 



A = mL 



(106) 



which measures the spatial variation in units of the 
Compton wave length, we obtain: 



d t +eE{t)d P3 -Miy s ,p) w«(x,p;t) 



-eE(t)A{x 3 ) 



a,, + 



P3 t 24A="P3 



d S 



w< '(p;t) . (107) 



This equation serves as the starting point for our analy- 
sis of the influence of a small spatial inhomogeneity. In 
order to estimate the influence of higher derivatives, we 



compare the terms occurring in the derivative expansion 
of the pseudo-differential operator D t in Eq. (|102p : 



and 



d P3 w {0) (p;t) with 



2fV^ 3 w(°)(p;t) 



(108) 



In fact, by means of this procedure we do not quantify the 
overall influence of the small spatial inhomogeneity A (#3) 
on the Schwinger effect; for this we would really have to 
solve the PDE system Eq. (|101[) . However, by means of 
the derivative expansion we might estimate a parameter 
region for which the higher derivatives do not play an 
important role such that we could restrict ourselves to 
the solution of the first-order PDE system: 

[dt + eE(t)d P3 ] w(x,p;t) = M(V s ,p)w(x,p;t) . (109) 

It is known from the analysis of the Schwinger effect in 
spatially homogeneous electric fields that the orthogo- 
nal momentum solely acts as an additional mass term 
and does not change the qualitative behavior. Thus, for 
simplicity, we will restrict ourselves in the following to 
p± = 0, such that we deal only with the following DHW 
functions: 



s(p3\t) 

V3(p3\t) 
til, 3(^3;*) 



W(P3) 
2^3 



x'Cfc;*) 



2p 3 



X Z (P3;t) , (110) 



w(p 3 )' 

x'ih-.t) + ^— x 2 (p3-,t) , (111) 



Uj(p 3 ) ' LJ{p 3 

2x 3 (P3;t) , 



(112) 



^,(0) = J_r_ 3? 3 3S (0) +2v (0) 



i n,3 



d 3 v (0) 

d 3 1 (0) 

U p 3 % l,3 



4 

7? 



4 
7% 



(0) , 2d, 2 (p 3 ) (0) 



b l,3 



(115) 



2p 3 u 1 2 (p 3 ) jQ) _ 2^ 2 (p 3 ) .,(0) 

e S i V 3 



3p 3 tn 



where we have dropped the arguments of the DHW func- 
tions for simplicity. 



According to Eq. (|113p . we compare dp. 



wW( 



with 



'ps**- '(Pa) 

^4<9? w^(ps) in the following; due to the fact that all 
the DHW functions show a rather similar behavior, we 

restrict ourselves to t^ 3(753) for simplicity. Fig. Q] clearly 
shows that the higher derivatives always become more 
important than the first derivative for large momenta p 3 . 
This might be understood in the following way: Due to 
the acceleration in the electric field, all length scales are 
ultimately probed and become important, even though 
the pair creation process happens on the length scale 
of 0(A). Note, however, that the point at which the 
higher derivatives are of the order of the first derivative 
depends on the electric field strength e = Eq/E c : For 
higher field strengths this point is already reached for 
lower momenta. 

In order to estimate the importance of the spatial inho- 
mogeneity for the pair creation process itself, we should 
have a closer look at the pair production rate in Fig. [5] 

For e = 0.2, the dominant contributions to the pair 
production rate stem from a region of kinetic momenta up 



where we have introduced the dimensionless phase-space 
kinetic momentum p 3 = — and the dimensionless energy 

variable w(p 3 ) = " (P3 - . 



B. Example 1: Constant electric field 

We explicitly calculated the time-independent coeffi- 
cients x i={1,2,3} (P) m App.E such that the DHW func- 
tions Eq. (|110p - (|112p for p j_ = and in terms of the 
dimensionless variables p 3 and £>(p 3 ) are easily calcula- 
ble. In the end, we want to compare the first with the 
third derivative of the source term in Eq. (| 1 T[) : 



-meA(x 3 ) 



"p3 + W>?®P3 



w(°>(p 3 ) 



(113) 



Therefore, we explicitly determine these derivatives ac- 
cording to Eq. (p4| . yielding: 




%s (0) = ^40) 



,(0) _ 



3 • 



(0) 



U P3 w 3 — n, 3 ' 



FIG. 1: Comparison of the leading-order derivative term 



0, 



s ti 3(^3) (dashed) with the next-to-leading order term 

(114) M^tii !^) (solid) for E(t) = E Q with e = E Q /E C = 0.2 
(upper) and e — 1 (lower). For higher momenta, the next-to- 
leading order term eventually exceeds the leading-order term 
as the electric field persistently accelerates the produced pairs. 
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emphasizes higher momenta such that higher derivatives 
always become more important than the first derivative. 



C. Example 2: Sauter-type electric field 

For the Sauter-type electric field we may perform a 
similar analysis like for the constant electric field. Again, 
we only consider p± = 0, however, there is a huge qual- 
itative difference since the DHW functions now depend 
on both the phase-space kinetic momentum p^, and the 
time variable u. Using the expressions given in App. [B] 
we are able to analytically calculate the DHW functions 
Eq. (IllOp - (|112p . Again, we consider the source term in 
Eq. (|107p . which now reads: 



4meu(l — u)A(x3) [dp 3 + 



i «3 
24A 2 ' U P3 



><°>(P8;«) (116) 



FIG. 2: Pair production rate for E(t) = Eo with e = Eo/E c — 
0.2 (upper) and e = 1 (lower). For e — 0.2, pair production 
occurs dominantly for momenta pz < 2 where next-to-leading 
order derivative terms remain small, cf. Fig. [l] By contrast, 
higher derivative terms are expected to take a quantitative 
influence on pair production for e = 1 and sizable spatial 
variation A. 



to ^3 ps ±2. In this regime, the first derivatives are still 
of the order of the third derivatives, as shown in Fig. Q] 
As a consequence, the effect of higher derivatives on the 
pair production process should be taken into account only 
for a spatial variation of the Conrpton wavelength A = 
0(1), whereas the effect of the higher derivatives becomes 
suppressed for larger variation scales A . 

For e = 1, the situation is slightly different: Again, the 
donrinant contributions to the pair production rate arise 
from kinetic momenta up to j>> w ±2. however, there 
are non-vanishing contributions for higher monrenta as 
well, as shown in Fig. [TJ Due to the fact that the higher- 
derivatives become more important than the first deriva- 
tives for large momenta, the scale of spatial variation A 
has to increase as well in order to suppress the higher 
derivatives. 

To conclude: Concerning the pair production pro- 
cess, there is a strong interplay between the electric field 
strength e and the scale of spatial variation A. This in- 
terplay between field strength and scale of spatial varia- 
tion affecting the quality of the derivative expansion has 
already been observed in earlier studies [38| . In order 
to keep the pair production process itself unaltered by 
higher derivatives, the scale of spatial variation A must 
get larger for higher electric field strengths e. However, 
even if the effect of the higher derivatives on the pair pro- 
duction process might be negligible, the effect on the final 
nronrentunr distribution might be large: this is due to the 
acceleration of the pairs in the electric field which finally 
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FIG. 3: Comparison of the leading-order derivative term 
dp a t\ l(p3; u) (dashed) with the next-to-leading order term 

l/24d$ 3 t ( °l(p 3 ;u) (solid) for e = E /E c = 0.2 and 7 = 1 for 
t — —t (upper), t — (middle) and t = r (lower). 
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In what follows, 
tive 9p 3 w^ \p3] u) 

f (0) 



'(pa;u). 



we compare the first deriva- 

with the third derivatives 

Again, we restrict ourselves to 



tii 3(^3; u). Nonetheless, there are three big differences 
in comparison to the constant electric field: First, we are 
not able to calculate the first and third derivative with 
respect to ^3 as easily as for the constant electric field, 
as shown in Eq. (|114[) - (|115p . since these derivatives 
now include parameter derivatives of the Gauss hyperge- 
ometric function. Nevertheless, a numerical calculation 
is rather simple. Second, the situation is not quasi-static 
anymore but it makes a difference at which moment of 
time the system is considered. Third, the field strength 
e is not the only relevant parameter but we also have to 
consider the dependence on the pulse length r via the 
Keldysh parameter 7 = l/(mer). In order to analyze 
the interplay between these quantities, we consider the 
system at three different instants of time: At t = — r 
(increase of field strength), t = (maximum of field 
strength) and t = t (decrease of field strength). 

Let us first concentrate on Fig. [3] with e = 0.2 and 
7 = 1, where we already observe some general features: 
As in the case of the constant electric field, the third 
derivatives show an out-of-phase behavior compared to 
the first derivative: At momentum values where the third 
derivatives show a local maximum, the first derivative 
shows a local minimum and vice versa. As a consequence, 
depending on the scale of spatial variation A, the influ- 
ence of the first derivative might be inverted due to the 
third derivative term. We also observe that the relative 
importance of the third derivatives in comparison to the 
first derivative becomes bigger for later times. This might 
be interpreted in the following way: At late times the cre- 
ated particles have had more time to be accelerated in the 
electric field and, as a consequence, have travelled over 
larger distances and were exposed to even large scale in- 
homogeneities. The effect of acceleration in the electric 
field might also be seen in the shift of the global maxi- 
mum of the third derivatives towards higher momentum 
values. If we want to suppress the higher derivatives 
in comparison to the first derivative, we should choose 
the scale of spatial variation to be at least of the order 
A>0(5). 

Next we switch to a longer pulse with 7 — > 0.2 as shown 
in Fig.[4l which corresponds to a pulse with the same elec- 
tric field strength but with a 5-times longer duration. We 
see that the behavior at early times t = — t is very sim- 
ilar, however, the situation changes drastically for later 
times: First, due to the longer pulse duration the created 
particles are accelerated to higher momenta. Addition- 
ally, we observe a strong enhancement of the oscillatory 
structure and a strong increase of the magnitude of the 
third derivatives. Especially at t = t, the third deriva- 
tives are orders of magnitude larger than the first deriva- 
tives. This means that the scale of spatial variation A 
has to be chosen much larger in order to suppress the 
influence of higher order derivatives. 

The change in the overall magnitude is the distinctive 



feature when we increase the field strength e — > 1 as 
shown in Fig. [5] The general behavior is rather similar 
as before, however, the source term becomes more im- 
portant in comparison to the left-hand side of Eq. (|107l) . 
As a consequence, we expect the overall effect of spatial 
inhomogeneities to be more important for strong electric 
fields than for weak electric fields. 

To summarize, Fig. [3]- Ogive us the following physical 
picture: 

(a) Concerning the time variable u, we see that the 
role of higher derivatives is more important at late 
times than at early times. This suggests that even if 
the influence of spatial inhomogeneities on the pair 
production process itself might not be very strong, 
the final momentum distribution can be strongly 





P, [m] 



t=T 
0.4 - 




FIG. 4: Comparison of the leading-order derivative term 
dp 3 til(p3;u) (dashed) with the next-to-leading order term 

l/24^ 3 ti^(ps;u) (solid) with e = E /E c = 0.2 and 7 = 0.2 
for t = — r (upper), t — (middle) and t — r (lower). Note: 
For t = r the third derivative is in fact 2 magnitudes larger 
than the first derivative but has been scaled by a factor of 
0.01. 
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FIG. 5: Comparison of the leading-order derivative term 
d$ 3 til(p3;u) (dashed) with the next-to-leading order term 

l/24dj ti°^(p3;n) (solid) with e = 1 and 7 = 0.2 for t = -t 
(upper), t — (middle) and t = r (lower). 



altered. 

(b) Concerning the field strength e, we conclude by 
comparing the magnitude of the source term for 
the different values of e that the influence of spa- 
tial inhomogeneities on the pair production process 
becomes more important for higher field strengths. 
We interpret this observation as arising from the 
fact that for weak fields the created particles are 
not accelerated that much and, as a consequence, 
they do only feel inhomogeneities on large length 
scales but not on shorter variation scales. 

(c) Concerning the Keldysh parameter 7 we note that 
small values of 7 correspond either to strong elec- 
tric field strengths or to longer pulse durations r. 
Therefore, particles are more substantially acceler- 
ated in the electric field and, as a consequence, spa- 
tial inhomogeneities might have a greater impor- 



tance, since the particles feel the inhomogeneities 
even on shorter scales of variation. 

As a consequence, we observe that there will gener- 
ally be a complex interplay between all the relevant pa- 
rameters for any type of space- and time-dependent elec- 
tric field E(x,t). Thus, it is difficult to predict a priori 
whether neglecting higher derivatives can be a good ap- 
proximation or not. It is clear that in the limit of a 
spatially homogeneous electric fields A — > 00 the leading- 
order derivative approximation Eq. (|10*j|) becomes exact 
since all higher derivatives are suppressed by a factor of 
A~ 2 ™. However, it is not clear a priori for which values of 
A the higher derivatives play a quantitatively important 
role and should be taken into account. Future investiga- 
tions on that problem should help to better understand 
the interplay between the different scales. 



IV. CONCLUSIONS & OUTLOOK 

We have investigated the Dirac-Heisenberg-Wigner 
(DHW) formalism for non-perturbative pair production 
in general electromagnetic fields. As a genuine real- 
time formalism, this approach provides for a comprehen- 
sive framework of addressing all aspects of pair produc- 
tion most notably, the non-equilibrium character of pair- 
production in a fully time- and space-resolved manner. 

We have shown that the DHW formalism includes 
quantum kinetic theory (QKT) which has so far been the 
most successful approach to describe the real-time evo- 
lution of pair production in the limit of time-dependent 
but spatially homogeneous electric fields. We conclude 
that the DHW formalism provides for the desired gener- 
alization of QKT to the case of arbitrarily general space- 
and time-dependent electromagnetic fields. For a given 
field, the solution of the DHW formalism is parameter- 
ized in spinor QED in four-dimensional spacetime by 16 
irreducible components of the Wigner function which en- 
code the phase-space distributions of physical quantities 
such as mass, charge and current densities as provided by 
the produced pairs. From the knowledge of these quan- 
tities, although they are no semi-positive definite proba- 
bility distributions, physical observables such as the pair 
distribution function in phase space can directly be in- 
ferred. 

Whereas the DHW formalism is completely general as 
far as the details of the external field are concerned, we 
have confined ourselves in the present work to an analysis 
of exactly solvable cases such as the constant electric field 
and the Sauter potential. Of course, such exactly solvable 
cases always provide for a controlled starting point for 
more general cases, in particular, they should also serve 
as a benchmark for future full numerical studies. More- 
over, we have used these cases in the present work to pro- 
vide for a first glance at the possible use and limitations 
of natural approximation schemes such as the derivative 
expansion. The leading-order of this approximation cor- 
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responds to a locally- constant field approximation, i.e., 
approximating the spatial dependence of the field locally 
by a constant field. 

The picture arising from this investigation is rather 
diverse: it has already been known, for instance, from 
worldline instanton studies [33], (34(, that the locally- 
constant field approximation underestimates the pair- 
production rate for fields varying in time, and overes- 
timates the rate for fields varying in space. Whereas the 
exact solution, e.g. for the Sauter potential in Sect. Ill C~2l 
reflects this fact, we observe that general statements 
about the potential quality of the derivative or locally- 
constant-field approximation cannot straightforwardly be 
made. For all concrete examples, we observe that the 
next-to-leading order derivative terms in fact exceed the 
leading-order terms either for higher momenta or for late 
times. Taken at face value, this seems to imply that the 
derivative expansion is always bound to fail as soon as 
the field exhibits spatial variations. However, the reason 
for this strong modification of the next-to-leading order 
terms lies in the fact that the persistent presence of ac- 
celerating field components, of course, exerts a strong in- 
fluence on high-momentum components which eventually 
resolve also small spatial variations of the field. There- 
fore, it is only natural to expect that the derivative ex- 
pansion should fail in the way it does for high momentum 
components and at late times. 

Nevertheless, our results also provide a guideline to a 
less strict view on the quality of the derivative expansion: 
phenomenologically, the most relevant quantity is the 
pair distribution function in momentum space. Quanti- 
tatively, the question needs to be addressed whether the 
next-to-leading order terms of the derivative expansion 
exert a strong influence on this distribution function. As 
the higher derivative terms become dominant for higher- 
momentum components, we conclude that the derivative 
expansion can still remain a reasonable approximation as 
long as the dominant pair distribution is peaked at lower 
momenta. Whether or not this is the case, depends not 
only on the scale of spatial variation of the field, but also 
on the overall field strength and also possible further time 
dependencies. As a rule of thumb, we observe that the 
dominant low-momentum components appear to remain 
little affected at significantly subcritical field strengths 
with spatial variations being substantially larger than the 
Compton wavelength. 

Beyond the technical question about the quality of the 
derivative expansion, the most interesting question is cer- 
tainly as to whether a phenomenologically relevant inter- 
play between characteristic signatures of pair production 
and space-time shaping of the external field exists. Based 
on the suprising observations that have already been 
made for simple time-dependencies [45|, |46|, we expect 
that this question can be answered in the affirmative. 
From our present studies of spatial variations and the 
dominant effect on high-momentum components, we con- 
clude that temporal and spatial pulse shaping can have 
substantial effects on the momentum structure of the pair 



distribution function. As to whether temporal and spa- 
tial pulse shaping can also optimize the total number of 
produced pairs certainly remains the most pressing ques- 
tion which we hope to address with the DHW formalism 
based on full numerical solutions in the near future. 
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Appendix A: Quantum Kinetic Theory (QKT) 

In this appendix, we give a brief derivation of the quan- 
tum kinetic equation describing Schwinger pair produc- 
tion in spatially homogeneous, time-dependent electric 
fields J2J, [25| . As in the main part of this paper, we 
adopt the temporal gauge Aq = 0. We choose the vec- 
tor potential A(t) = A{t)e^ such that the Dirac equation 
reads: 

(ij°d t + ij- V s -ieA(t)e 3 -m]f(x,<) = 0. (Al) 

Due to spatial homogeneity, we decompose the spinor 
field into its Fourier modes according to Eq. (|29[) . such 

that the Dirac equation for the mode function ip(q,t) 
reads: 



(ij°d t - 7 • n(q, t) - m) tp(q, t) = 



(A2) 



Again note that T?(q,t) = q — eA(t)e 3 denotes the time- 
dependent kinetic momentum on the trajectory whereas 
q denotes the canonical momentum. In order to solve 
this equation we apply the ansatz: 

$(q, t) = (z 7 °9 t - 7 ■ n(q, t) + m) 4>{q, t) , (A3) 

such that the spinor-valued function cf>(q, t) obeys the 
equation: 



(d 2 t + £ 2 (<f, t) + leEit)^^) <t>(q, t) = , 



(A4) 



with uj(q,t) defined as before. It is convenient to expand 
</>(q,t) in a basis consisting of the eigenvectors of 7°7 3 , 
such that: 

$(q,t)=J2R 8 g s (q,t) with 7 7 3 ^ = Ai? s . (A5) 



There are two eigenvectors R s =i,2 with A = +1 and two 
eigenvectors R s =3,4 with A = — 1. Inserting this ansatz 
into Eq. (|A4[) . each g s (q, t) obeys the equation of a time- 
dependent oscillator: 



(df + w 2 (q, t) + ieE(t)X) g s (q, t) = 



(A6) 
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which are in general not exactly solvable; exceptions are 
the constant electric field E(t) = E and the Sauter-type 
electric field E(t) = E scch 2 (t/r) (see Sec. IITCJ) . Each of 
them is a second-order differential equation and possesses 
as such two linearly independent solutions g s (q, t) . Due 
to the fact that Eq. (|A4|) gives four equations but the 
ansatz Eq. (|A5|) allows for eight solutions, there is a re- 
dundancy which is removed by choosing only one set of 
eigenvectors, either s = {1,2} ors = {3,4}. Due to the 
absence of magnetic fields, we impose the same initial 
conditions for both spin states, such that: 



;(±) 



:(±) 



gi ( *>(q,t)=g 2 [ *>(q,t)=g^ 



(Q,t) 



(A7) 



Consequently, we canonically quantize ip(q,t) accord- 
ing to Eq. (|30p by introducing anti-commuting cre- 
ation/annihilation operators as well as four-spinors, with: 

u s (q, t) = {ifdt - 7 • n(q, t) + m) gM(q, t)R s , 

(A8) 

v s (-q, t) = (i^dt - 7 • Tf(q, t) + m) g { -\q, t)R s . 

(A9) 

Due to the fact that we work in the Heisenberg picture, 
the creation/annihilation operators are time-dependent 
in general, however, due to the choice Eq. (|A8[) - (|A9|) . 
the whole time-dependence can be absorbed into the four- 
spinors. It can be shown, that in the case of vanishing 
electric fields, the properly normalized vacuum solutions 
are given by: 



i vac 



(?,*) 



D Ti^(q)t 



v / 2uj(q)(uj(q)Tq3) ' 



(A10) 



with u)(q) 



\fm 2 



It is important to note that 



a particle/antiparticle interpretation of the field quanta 
is only possible in the case of such plane-wave solu- 
tions. However, as soon as electric fields are present, the 
mode functions g( ' (<f, t) are no plane waves anymore 
and an interpretation in terms of particles/antiparticles 
is not straightforward. It is a further consequence of the 
presence of electric fields that the Hamiltonian operator 
achieves off-diagonal elements which account for parti- 
cle/antiparticle creation/annihilation. 

The Hamiltonian operator might be diagonalized by 
performing a unitary non-equivalent change of basis to 
a quasi-particle representation via a time-dependent Bo- 
goliubov transformation: 

a 3 (q,t) = a(q,t)a s (q) - P*(q,t)bl(-q*) , (All) 

bt(-q,t) = mt)a s (q)+a*(q,t)bt(-q) , (A12) 

with the creation/annihilation operators becoming time- 
dependent but still fulfilling the equal-time anticommu- 
tation relations. In order to be a canonical transforma- 
tion, the Bogoliubov coefficients a(q,t) and f3(q,t) have 
to fulfill: 



In pure vacuum when no electric fields are present, 
the two different operator bases coincide such that 
a vac (q,t) = 1 and f3 vac (q,t) = 0. Note, that this relation 
also holds in the presence of electric fields at asymptotic 
times t — >• — oo. Within this so-called adiabatic basis, the 
Fourier modes ip{q,t) read: 

$(q, t) = ^U s (q, t)a s (q, t) + V s {-q, t)b\ (-£ t) . 

(A14) 
The adiabatic four-spinors are chosen in close analogy to 
Eq. (|A8[) - (|A9[) such that they coincide with the vacuum 
solutions in the case of vanishing electric fields: 

U s (q, t) = ( 7 w(<f, t) - 7 • *?(£ t) + m) G« (q, t)R s , 

(A15) 

V s (-q, t) = (-7°2(g, t) - 7 • 7r(q, t) + m) &^ (q, t)R s , 

(A16) 

with the adiabatic mode function G^ (q, t) given by: 

e =Fi0(g,*o,t) 



G^\q,t) = 



y/2w(q,t)(Lj(q,t)TMQ3,t)) ' 
and the dynamical phase Q(q,tQ,t) being defined as 



(A17) 



Q(q,to,t)= / u{q,t')dt' 



(A18) 



to 



The lower bound to is not determined since it only fixes 
an arbitrary phase at a given instant of time. Note that 
an interpretation in terms of particles/antiparticles is 
only straightforward at asymptotic times when the ex- 
ternal electric field vanishes and the solutions Eq. (JA17J) 
behave like plane waves. 

In order to define the single-particle momentum dis- 
tribution function f(q,t), we assume that we start with 
vacuum initial conditions at t — > — oo: 

(0| ai(q)a s (q) |0) = (0| b\(q)b s (q) |0) = . (A19) 

We then define f(q, t) as the instantaneous quasi-particle 
number density for a given canonical momentum q. Due 
to the absence of magnetic fields, we take the sum over 
both spin states, such that: 



f{q,t)= lim Y, 



V— >oo 



ai(q,t)a s (q,t) 
V 



2\(3(q,t)\ 2 ■ 



Hq,t)\ 2 + |/3(g, t)| : 



1 



(A13) 



(A20) 
with V being the (infinite) configuration space volume. 
As a consequence, the knowledge of f3(q, t) allows for the 
calculation of f(q,t). In fact, the different representa- 
tions of ip{q,t) Eq. (|30|) and Eq. (|A14|) translate into an 
expression for the Bogoliubov coefficients: 

a(q,t) = ie ± &-Hq,t) [d t - iu{q,t)]f+> >(q,t) , (A21) 
mt) = - t e ± & + \q,t)ld t + iu(q,t)]g(+Hq,t) , (A22) 
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such that their time derivatives form an ODE system: 

(A23) 



f t a(q,t) = lQ(q,t)/3(q,t)e 2 ^^ 



j t mt) = -\Q{q,t)a(q,t)e-^^ , (A24) 

with Q(q,t) defined as in Eq. (|4"T|). Introducing C(q,t) = 
2a(q, t)(3*(q, t), this ODE system might be rewritten as: 



Jt^ l) 



Q(q,t) l-f(q,t) 



-2iB(q,t -t) 



^M*)= -<9(?>*)Re CCfctJe 2 * 6 ®*"*) 



(A25) 
(A26) 



Formally integrating the first equation from a time of 
pure vacuum i vac — > — oo to t, yields the Vlasov equation 
for the single-particle momentum distribution function 
f(q,t) in its integro-differential form: 



<// 



mt)=Q(q,t)ff dt'Q(q,t>) 



x[l-f(q,t')} cos [2S(q,t,t)], (A27) 



with /(<f, i V ac) = 0. It is possible to rewrite this integro- 
differential equation in terms of an equivalent ODE sys- 
tem by introducing auxiliary functions p(q, t) and a(q, t): 



dt 

d_ 

~di 
d_ 

dt 



f(q,t) = Q{q,t)p(q,t) > 



(A28) 



p{q,t) = Q(q,t)[l- f(q,t)]-2u(q,t)Z(q,t) ,(A29) 



a(q,t) = 2u(q,t) p(q,t) , 



(A30) 



with appropriate vacuum initial conditions p(q, t v 
a{q,t vac ) = 0. 



Appendix B: DHW functions for exactly solvable 
electric fields 



constant electric field and the Sauter-type electric field. 
The DHW functions are obtained as follows: We already 
derived the analytic expressions for the single-particle 
momentum distribution function f(q, t) in Sec. Ill Cl As a 
consequence, according to Eq. (|4"8"|) - ([50)) we are able to 
calculate the non- vanishing coefficients X* = ^ 2 ' 3 HSj*)i c f- 
Eq. (|46|) . as well. Finally, we obtain the DHW functions 
according to Eq. (|5Tj) - (|55|) after performing the phase 
space projection Eq. (|4"0)) . 



1. Constant electric field 



In order to simplify the expression for the single- 
particle momentum distribution function f(u) derived in 
Eq. (|T2")) , we introduce the following abbreviations: 



di(u) 

d2(u) 
d 3 (u) 

d±(u) 
which fulfill: 



\ D -i+i V /2 {-ue l i 

i .• IT \ I 2 



D 



II) 



/2 {-ue °\ 



(Bl) 
(B2) 



■*■ x 



e t *D_ 1 _ iv/2 (-ue~ 

D l7l/2 {-ue-^)+c.c. (B3) 
e~ l ^D_ 1 _ lv/2 (-ue'^) x 

D lv/2 (-ue- l2 i) + c.c. (B4) 



d u di(u) 
d u d 2 (u) 
d u d 3 (u) 



d^u) 
-fd 4 (u) 
-udi(u) 



(B5) 
(B6) 
(B7) 



In this appendix we give the analytic expressions for 
the DHW functions for the exactly solvable cases of the 



J 



We may then express f(u) in terms of di(u), such that 



A' 



:j={l,2,3} 



(u) are given by: 



^-^IHra 



di(u)+ 1 + 



VV 



d 2 (u) 



V^T 



,d 3 (u) , (B8) 



X 2 (u) 
X 3 (u) 



V 1_ 

2 V^+t 



= e 4 



? {-r)d 1 (u) + 2d 2 (u)+ud 3 (u)\ , (B9) 

(2* + «»)>/» V2(2,TW^ e ""^" l(U) ~ 2 " 2H + (2?? + U2)3/2d>) } ' (B10) 



2 '/ 



r 



with x l { u ) = 1 — /( u )- I n order to obtain the DHW p + eA(t). Due to the linear relation between 93 and t, 
functions we perform the variable transformation q —> 



this phase-space projection is trivial and reads: 



P3 



u -> \/7P3 and di(p) = d t 



(Bll) 



where we introduced the dimensionless phase-space ki- 
netic momentum p 3 = —■ Note that the di(p) implic- 
itly depend on the orthogonal kinetic momentum p± by 



J 
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means of r\ — e\/eEo with e\ = m 2 + p\. Obviously, 
di(p) do not depend on the time variable t but only on 
the kinetic momentum p such that the Schwinger effect 
in a constant electric field might be regarded as a quasi- 
static problem. The phase-space coefficients p^f 1 - 2 * 3 } (p) 
which allow for the calculation of the DHW functions 
Eq. (J5TJ) - ([HHD thus read: 



x\p) 
x 2 (p) 
x 3 (p) 



1 - e 4 " E o 



2eE 



- -^r 1 d x {p) 



u(p) 



P3, 

uj(p) 



d 2 (p) - 



V^Eo^ip) 



d 3 (p) 



eE ej_ 
2uj 3 (p) 



.-±-d 1 (p) + 2d 2 (p) 
eE 



eEi 



■P3d 3 (p) 



1 + -e 4eE o 



-^-d 1 {p)-2d 2 {p)+(^ r 
eE \ eE 



3/2 



w 3 (p)d 4 (p) 



(B12) 
(B13) 
(B14) 



2. Sauter-type electric field 

We start from the expression for the single-particle mo 
mentum distribution function /(<?3,u) given in Eq 
and introduce the following abbreviations: 



hi(q 3 ,u) = F(a,b,c;u) 



(B15) 
(B16) 



M<73,u) = fF(l + a,l + b,l + c;u) 

h 3 {q 3 ,u) = -ifF(l + a,l + b,l + c;u)x 

F(a*,b*,c*;u) + c.c. (B17) 



and 

Qi(Q3:U) = p(q 3 ,u) - (1 -u)Q(q 3 ,0) -uw(g 3 ,l)] 2 (B18) 
&(«3,«) = 4 7 2 eV(l- U ) 2 (B19) 

Q 3 (q 3 ,u) = 27e«(l-«) x 

[w(fla,«) - (1 - u)w(«s, 0) - «£)(&, 1)] (B20) 

such that f(q 3 ,u) can be written as: 



f(q 3 ,u)=N f (q 3 ) 1 + 



^3(93, -») 

W(?3,«) 



^ 0i(?3,«) /li (©,«)■ 

(B21) 



Again note that {a, 6, c}, which have been defined in Eq. (|85[) . only depend on 53 but not on u. Taking into account 
the general derivation formula for the Gauss hypergeometric function Eq. (|96p . we can explicitly calculate the first and 
second derivative of Eq. (|B21|) . After calculating d u f(q 3 ,u) and d 2 f(q 3 ,u), we are able to determine the coefficients 



X 



j={l,2,3} 



(q 3 ,u) according to: 

x 1 ^,^) = 1 - f(q 3 ,u) 



X 2 (q3,u) 



jQ 2 (q 3 ,u) 
2V1+K 2 



d u f{q3,u) 






7S (g3.")^3(g3,") 
1 + K 2 



^7(g3^)- 7 Igjl\']' ) ^/(g3,u)} 



(B22) 
(B23) 

(B24) 



r 



In order to obtain the coefficients in phase-space, we variable transformation: 
have to perform the variable transformation q — > p + 
eA(t), which reads: 



^3(93, u) -> p 3 , 



q 3 -> p 3 



(B25) 



Q(q 3 ,u) -> w(p 3 )=y/l + K 2 +P3 



(B26) 
(B27) 



As a consequence, the quantities Q(q 3 ,u) and 7T3(<?3,u) 



On the other hand, any function of the canonical mo- 
mentum q 3 only, e.g. tD(g 3 ,0) or £((73, 1), acquires a de- 



only depend on p 3 but not on u after performing this pendence on both the phase-space kinetic momentum p 3 
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and the time variable u: 



and 



TrsfeO) = q 3 - i ->&- *f , 



^3(93, 1) = 93 + ^ 



P3 



2u 

7 
2u-2 

7 



(B28) 
(B29) 



2( 9 3,0)^ \Jl + ^+(p 3 -^) 2 , (B30) 

Q(83,l)-^ A/l + ^ + fe-^V (B31) 



Therefore, whereas the functions hi(qs,u) depend on u 
solely through the last argument of the Gauss hypergeo- 
metric function, the transformed functions hi(p3; u) have 
a twofold u dependence: On the one hand, there is still 
the u dependence due to the last argument. On the other 
hand, due to the fact that the parameters Eq. (|85j) were 
function of 173 only, they will depend on both ^3 and u 
after the transformation to phase space. 
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